arXiv:nucl-th/9808022v3 24 Nov 1998 


Removing Pions from Two-Nucleon Effective Field Theory 


James V. Steele and R. J. Furnstahl 
Department of Physics 

The Ohio State University, Columbus, OH 43210 
(August, 1998) 


Abstract 


Two-nucleon effective field theory (EFT) beyond momenta of order the pion 
mass is studied for both cutoff regularization and dimensional regularization 
with power divergence subtraction (PDS). Models with two mass scales il¬ 
lustrate how effects of long-distance pion physics must be removed from the 
coefficients that encode short-distance physics. The modified effective range 
expansion shows that treating pions perturbatively, as in the PDS power 
counting, limits the radius of convergence of the EFT. Predictions from both 
regularization schemes with one-pion contributions are compared to data. The 
breakdown of the effective field theory occurs for momenta of order 300 MeV, 
even using the modified effective range expansion. This behavior is shown to 
be consistent with that expected from two-pion contributions. 


I. INTRODUCTION 


Nucleon-nucleon effective field theory (NN EFT) without pions, which applies at low 
energies, is now well understood. Predictions are made by constructing the most general 
Lagrangian with nucleon helds only and hxing the coefficients by matching calculated and 
experimental observables in a momentum expansion. A key feature of NN EFT is the need 
for a nonperturbative treatment, as implied by the existence of shallow bound states such 
as the deuteron. Heavy propagating nucleons lead to an enhancement of loop graphs [|^, 
which destroys the systematic power counting characteristic of perturbative EFT’s [0. 

Weinberg’s solution was to sum certain diagrams to all orders in the loop expansion via 
the Schrodinger equation |]l|,|^. Recent analyses hnd this summation to be problematic for 
systems with large scattering lengths [Q-Q , which has led to alternative prescriptions within 
various regularization schemes 0-^. The use of a given regularization scheme has implica¬ 
tions on the renormalization scheme chosen because in practice calculations are truncated. 
The subsequent calculations of low-energy two-nucleon amplitudes have been shown to be 
systematic and regularization scheme independent |]IO i pT| , p|, reproducing the effective range 
expansion of nonrelativistic potential scattering 0. The unambiguous predictions of addi¬ 
tional observables are in good agreement with data, as seen, for example, in the triplet state 
of threshold nd scattering |^, which is dominated by two-body scattering. 

The characteristic breakdown scale, or radius of convergence, for a momentum expansion 
of the low-energy NN R-matrix is dictated by the pion. In conventional effective range theory, 
this breakdown is manifested by a branch point in pcotS at p = m,r/2 from the one-pion- 
exchange Yukawa potential (see Sec. II). More generally, the physical picture from effective 
held theory anticipates a breakdown for momenta of order p ~ as details comparable to 
the Compton wavelength of the pion are resolved. Predictions for external momenta above 
this point require taking the pion explicitly into account. In this paper, we address some of 
the open questions about power counting with pions and the radius of convergence of the 
NN EFT. 

In Ref. 0, a power counting scheme was introduced, in which all pion contributions to 
the scattering amplitude enter perturbatively. Calculations were simplihed by the use of 
dimensional regularization with power divergence subtraction (PDS). In contrast to Wein¬ 
berg’s scheme, this new approach sums only the leading contact interaction between nucleons 
to all orders. 

Forcing a summation of higher-order contact interactions can be devastating in some 
regularization schemes — such as dimensional regularization with minimal subtraction — 
leading to undesirable systematics for a large scattering length HJTT|,0. PDS was intro¬ 


duced to dehne a consistent power counting scheme that avoids this forced summation. 
However, with a regularization that is solved numerically, such as cutoff regularization, it 
is simpler in practice to sum up these higher-order terms. We have shown previously that 
regularization schemes exhibiting naturalness in the bare parameters have the proper EFT 
behavior even when higher-order terms are summed O- In other words, as long as the 
cutoff is suitably chosen so that the theory is natural (bare parameters of order unity), the 
summation of higher-order terms does not affect the power counting implying that 
cutoff regularization is an equally valid scheme for calculations. 

Properly adding the pion should, in principle, increase the range in momentum for which 
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the EFT is predictive to the mass of particles not explicitly taken into account in the 
Lagrangian. Optimistic estimates in the literature anticipate a breakdown on the order of the 
p mass or p ~ 800 MeV . However, both dimensional 0 and cutoff [Q regularization 

calculations of NN scattering with pions imply a breakdown closer to p ~ 300 MeV. The 
subsequent application of EFT techniques to nuclear matter could be problematic due to 
the close proximity of the Fermi momentum (280 MeV) to the breakdown point. 

The observed breakdown of NN EFT could have various sources. Explicit inclusion of 
correlated two-pion exchange, pion production, or the A(1232) may be necessary to increase 
the radius of convergence beyond several hundred MeV. The breakdown scale could even be 
mp/2. However, we should hrst establish that the EFT is indeed limited by physics and not 
by how it is implemented. In the cutoff approach, simply adding one-pion exchange to the 
effective potential might not completely remove its influence on the radius of convergence. 
Qualitatively different reasons may limit the PDS approach. A renormalization group anal¬ 
ysis showed that PDS power counting should break down at a new scale Ann ~ 300 MeV, 
but there has been no systematic error analysis for PDS with perturbative pions. 

In this paper, we address these implementation issues. Our approach generalizes the 
procedure of Ref. by requiring the consistent removal of long-distance (pion) physics from 
the coefficients that encode short-distance physics. In Sec. H, we illustrate a technique to do 
this, called the modihed effective range expansion [^1, which treats pions nonperturbatively. 
We demonstrate its effectiveness and relevance within the context of two toy models using 
cutoff regularization. The modihed effective range expansion is used in Sec. HI to show that 
treating pions perturbatively within the PDS scheme limits the radius of convergence of the 
EFT. In Sec. IV, we compare NN scattering data to results from applying these techniques 
to both regularization schemes. The breakdown around p ~ 300 MeV persists, but this scale 
is shown to be consistent with that expected from two-pion contributions. Section V is a 
summary of our results. 


II. EXAMINING THE NEED FOR THE 
MODIFIED EFFECTIVE RANGE EXPANSION 


The effective range (ER) expansion for nonrelativistic NN scattering 
using analyticity arguments to relate the momentum times the cotangent 
pcot5 to an expansion in energy E = p‘^/M, with M the nucleon mass,[| 


m is derived by 


of the phase shift 


pcot6 



oo 2n 

I 2 P 

n=0 


1 

2 


( 1 ) 


Here the expansion is written in terms of the lightest mass scale A that characterizes the 
underlying dynamics |^. For example, a Yukawa or exponential potential with mass A 
limits the region of analyticity and hence validity of Eq. (0 to p < A/2 m- Although the 


^We restrict the discussion throughout this paper to interactions in the center of mass between 
two nucleons with a reduced mass M/2, M = 940 MeV, and neglect electromagnetic effects. We 
also take rriT^ = 140 MeV, irip = 770 MeV, /tt = 92 MeV, and = 0.075 below. 
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scattering length can take on any value, we expect the ranges to be of order 1/A 
this is empirically shown below. Special considerations are needed for the behavior of these 
ranges when the underlying potential is strong, as discussed in Sec. IV. 

In our earlier work |^, we used a short-range delta-shell potential at r = l/rup in the 
radial Schrodinger equation as a clean laboratory with which to compare various regulariza¬ 
tion schemes. As each new order in the momentum expansion of the effective potential was 
matched to data, we verihed the improvement of the error 


A(p cot S) = \p cot 5efr — P cot ^tmel (2) 

by powers of p^. Contributions from different orders in the expansion become comparable 
at what is known as the radius of convergence of the effective held theory. We expect this 
to occur for momenta comparable with the underlying scale A. We demonstrated this by 
looking at error plots [Q of A(pcot(5) as a function of momentum. Both cutoff and PDS 
regularizations were found to give systematic behavior independent of the details of the 
underlying theory. 


A. Delta-Shell and Exponential Model 


We now add a long-range part to the underlying toy theory to model the inclusion of 
the pion in the analysis of NN scattering. We hrst study the addition of an exponential well 
II with mass m^r. 


777 

V{r) = - ''^0 e” 


( 3 ) 


The exponential and delta-shell model the long- and short-range potentials respectively. 
Our ability to solve Eq. (H) analytically allows us to validate the numerical procedures used 
below in more realistic cases. The exponential results in singularities in p cot b starting at 
p = —imjr/2 analogous to the branch points from a Yukawa [jT^. We use a “chiral-symmetry” 
coupling Vo = m^/M. The delta-shell coupling will be taken to be 5 ^ = 1.01 for comparison 
to our earlier paper ||rT| . 

We determine the phase shift for V{r) by calculating the A-wave dost function JT(p), 
which plays a key role in the modihed effective range expansion [Q. This is accomplished 
by solving the Schrodinger equation for the dost solution f{p,r), 


dr'^ 


-MV{r) 


f{p,r) = 0 . 


( 4 ) 


The dost solution satishes the asymptotic condition f{p, r —> cx)) —>• and has the property 

f{p, 0) = J^{p). The exact dost function for an exponential potential [1^ can be generalized 
for the potential in Eq. @ and written in terms of Bessel and gamma functions. 


np) = 


f)T(i 




J-u{yo) + -^^ — J-y{y) {Jy{y)J-y{yo) - J-y{y)Ju{yo)) 


sm TTZ/ TUt, 


( 5 ) 


with z/ = 2ip/mT,, = 2-\/MI^/m7r, and y = This satishes the well known 

properties of the dost function, namely = bF*[—p*) and for V ^ 0 or p —>• cx), JF(p) ^ 1. 
We can extract the exact phase shift dtrue by using JE(p) = 
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Expressing the Bessel functions in Eq. (|]) as a series, the Jost function can be expanded 
in powers of l/rUp, 

with Q = {m^,p} and 


^ £ «(1 - 0 


(fc-z/) 


(7) 


Equation (|^ shows explicitly the singularities starting at p = —imT^/2 that limit the region 
of analyticity of the effective range expansion. These singularities always enter multiplied 
by the strength of the potential yo and have less impact on the radius of convergence of the 
EFT as yo decreases (see below for examples). 

Based on the above discussion, using a purely short-range effective potential to describe 
V{r) should result in a breakdown at p ~ m^/2. For S'-waves, the effective potential is 


Kfr(r) = Co6^{r) - ^^2 + ... , 


( 8 ) 


which in momentum space becomes 


Ves{p,p) =Co + + . . . . (9) 

It can be hxed to a given order in p by matching the C 2 nS to the ER parameters, Eq. (§• 
We do this by guessing an initial value for each of the C 2 nS and calculating pcot(5eff. Then 
a £t of the difference 


A{pcot6) = a + (3 ^ + ... , 


( 10 ) 


to a polynomial in p^/A^ results in the constants a,P ,..., which are minimized up to the 
order of I4fr(^") with respect to variations in the C 2 nS [M. Weighting the £t by both the 


expected theoretical error and experimental uncertainty ensures that the EFT constants are 
properly obtained. The singular nature of I4fr(?") requires regularization. We concentrate 
on using cutoff regularization with a cutoff Ac in this section, as detailed in Ref. |^. We 
focus on PDS in the next section. 

The cutoff Ac is roughly adjusted to minimize the error, which identifies the resolution 
scale of the underlying short-distance physics ^j. The results for using a short-distance 
effective potential to describe Eq. (||) are only weakly dependent on Ac, but are best for 
Ac ~ 70 MeV, which are shown as the dashed lines in the left plot of Fig. 2. We plot the 
error A(pcot 6) as a function of momentum on a log-log plot. The lines with increasing slope 
refer to one, two, and three constants hxed in Ws respectively. The radius of convergence can 
be read from the graph to be about 90 MeV. This breakdown point is directly connected to 
the singularity in p cot 6 at m 7 r/ 2 , as can be conhrmed by varying the mass in the exponential 
ofEq. m. 
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FIG. 1. The model of a delta-shell plus exponential potential, Eq. (^), with low energy 
constants ht by using the conventional effective range expansion (left) and modified effective range 
expansion (right). The radius of convergence on the left is about 90MeV regardless of whether 
the exponential potential is explicitly included (solid) or not included (dashed) in the effective 
potential. Using the modified effective range expansion increases the radius of convergence to 
about 1 GeV. 


It might be expected that by including the long-range behavior of the exponential po¬ 
tential —explicitly in the effective potential Eq. (^, we can increase the radius 
of convergence beyond m7r/2. Although the values of the C' 2 n’s from matching to the ER 
parameters in this case change slightly, the subsequent error plots are almost unchanged, as 
shown by the solid lines in the left plot of Fig. |l|. Naively taking the long-range physics into 
account in the cutoff effective field theory is simply not enough to yield an improved result. A 
nonperturbative treatment of long-distance physics with a cutoff leads to contamination in 
the short-distance coefficients. This can be hxed by appropriate counterterms as discussed 
in Sec. III.C. But hrst, we show how to £x this problem using the modihed ER expansion. 


B. The Modified Effective Range Expansion 


Properly accounting for the long-distance physics (the pion) is equivalent to removing 
its effects from the ER expansion. This should result in a new function with a region of 
analyticity larger than pcotd and lead to an improved radius of convergence for the EFT. 
After separating a general potential into a long- and short-range part V = Vl + Vs, such 
a function can be constructed and expressed in terms of the full phase shift 6 and the dost 
function J^l{p), dost solution fLip,^), and phase shift 5l for the long-range potential only 

II: 


K{p) 


pcot(d-dL) r/L(T,0) 

\:Fl{p)? ^ [ -Fl(p) 


1 1 “ ^ p2n 

—— + -p y r. 


de 


n=0 


/^2r. 


( 11 ) 
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The prime denotes differentiation with respect to r. Equation (pT]), which is called the 
modihed effective range expansion, has been most commonly used to remove Coulomb effects 
from the analysis of NN scattering; it can also be generalized to higher angular momenta 
I l]!^. The combination K{p) removes the effect of the long-range potential Vl, increasing 
the relevant A to the order of the short-range potential. Thus we expect that hn ~ 1/A. For 
the toy theory of Eq. (|^), this would mean removing the pion singularity and moving the 
radius of convergence up to rup, the scale of the delta-shell. Note that when the long-range 
potential is zero, 6l = 0, /l(p, r) = and the conventional effective range expansion is 
recovered. 

We should emphasize that the reason we introduce the modihed ER expansion is to 
ensure that the addition of the long-range potential to our EFT does not contaminate the 
process of htting the low energy constants. Once the constants are ht properly, we can return 
to the evaluation of the familiar observable p cot 6 to assess the improvements made by the 
effective held theory. When the phase shift is close to zero, pcot 6 becomes abnormally large 
and it can be advantageous to look at error plots of K{p) instead. 

We now apply the modihed ER procedure to the toy model of Eq. (j|). We again include 
the exponential —in the EFT, but now match the low-energy constants C 2 n by 
enforcing AK{p) = 0 to the order in we are working, similar to what was done for 
A (p cot (5) in Eq. O- The exact dost function Eq. (H) with p = 0 is the long-distance 
dost function J^l{p) used in Eq. ([TT|). Since both the EFT and the true theory have the 
same long-distance behavior (by construction), the term in the modihed ER function that 
depends on //(p, 0) cancels in the diherence. 

The results of the modihed ER ht are shown in the right plot of Fig. The optimal 
(natural) value of A^ is now about 800 MeV and the radius of convergence is seen to be about 
1 GeV. Furthermore, the EFT constants that are matched using the modihed ER procedure 
agree exactly with the values from the delta-shell potential alone |^. This verihes that 
the long-distance physics is successfully removed. For more general potentials, with long- 
and short-distance scales of m-jr and rrip respectively, it can be shown that the modihed ER 
expansion Eq. (|Tl|) is the same as the pure short-distance ER expansion up to terms of 
O/rn^p) []^. These corrections are actually zero for a delta-shell potential. 

When comparing to data, we expect to have some error in the long-range potential used 
in our ehective theory compared to the underlying theory. We can estimate the ehect this 
will have on our results by purposely detuning the value of Vq used in the long-range part 
of our ehective potential. The ht to three constants is sensitive enough to fail if the error is 
larger than a fraction of a percent. The ht to one constant gives consistent results regardless 
of the error. The comparison of the results for a ht to two constants is shown in Fig. for 
errors of 5, 10, and 50 percent. Although the results are not as clean as in the error-free 
case, we can read oh a radius of convergence by taking a straight line from the low energy 
behavior. This shows that the radius of convergence inferred from the plot slowly reduces as 
the error in Vq increases. As long as the discrepancy in the long-range part of the ehective 
held theory and that from data do not diher by more than 5%, the results for one and two 
constants are fairly stable. 

Using the modihed ER expansion is a guaranteed way to properly incorporate the long¬ 
distance physics in the EFT. However, at least some improvement in the radius of con¬ 
vergence was found in the case of NN scattering without this technique PJ^. We use a 
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FIG. 2. The dependence of the two-constant fit to Eq. (^) on the percentage error in Vq. 

second model in the next snbsection to fnrther investigate when the modihed ER expansion 
is needed. 


C. Two-Yukawa Model 

A more representative toy model of nuclear forces consists of two Yukawas of masses 
and nip. We require 


V{r) 


(y.r 


0^7, 


( 12 ) 


but for small r the results are practically independent of whether a cut-off Yukawa or a 
hard core repulsion are used to regulate the potential at the origin. For simplicity, we take 
Eq. (^) as the exact potential and merely exclude the origin. The couplings are initially 
taken to be = 0.075, as given by one-pion exchange in NN scattering, and ap = 1.05, 
which gives the large scattering length found in that channel, = —25fm. Equation (p!2D 
models the short-distance physics with a single scale. Nuclear phenomenology suggests that 
the short-range Yukawa actually should be replaced by the combination of attractive and 
repulsive Yukawas resulting from the exchange of scalar and vector mesons [|T^. We will 
come back to this issue in Sec. IV. 

Since the long-distance physics is dictated by the pion, the effective Lagrangian that 
describes this model problem is identical to that for NN scattering lHIJ 


= Alt 


idt + — -■ VItt ■ r) 

2M 2U 
1 


N 


-CoiN^Nf + -Ga [(YtVV)2 + ((VYt)^)^] + ... , 


(13) 
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FIG. 3. The error in pcot5 (left) and the phase shift in degrees (right) for the two-Yukawa 
toy model. The solid curves on the left show errors from the modified ER fit to one, two, or three 
constants, which are independent of the strength of Vl- Without using the modified ER procedure, 
similar results are obtained for one and two constants with a weak The dashed curves on the 

left show the errors as well as the radius of convergence worsen if is made a factor of three 
stronger. 

with the ellipses representing terms with more than one pion field as well as higher-derivative 
terms. Neglecting retardation effects since they are absent from the underlying toy model 
by construction, Eq. ([T^) can be expressed in terms of a potential (for the ^S'o channel) 

V,s{r) = 5\r) - 1^2 {yH\r) + + ... , (14) 

with = g\rril/167if^, showing that the long-distance physics of the pion from Eq. ([I^) is 
included exactly. The full effective potential, including the pion Yukawa, is regularized at 
short distances by a cutoff Ac [§]. 

The dost functions for both the full phase shift and the long-distance part required in 
Eq. (|T^) are calculated numerically, as outlined in Appendix A. Note that when calculating 
the modihed ER expansion K{p), we still impose a cutoff on the potential Eq. ([T^), as would 
be done if the underlying potential were not exactly known. 

The results from using the modified ER expansion to £t one, two, and three constants 
are shown by the solid lines in the left plot of Fig. for the optimal value Ac = 600 MeV. 
The phase shifts themselves are shown in the right plot. While the improvement in the 
phase shifts is manifest, it is clear that the error plots are needed to make a quantitative 
assessment of the convergence. The plots show the radius of convergence is between 500 
and 600 MeV, which is comparable to the range of analyticity mp/2 of the modihed ER 
expansion. For comparison, we also used a delta-shell at r = l/rrip as the short-distance 
potential (see Appendix A). This gives a radius of convergence of about 1 GeV, which is in 
agreement with the exponential potential model in the previous subsection, and shows the 
different impact from branch point and pole singularities in p cot 6. 
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TABLE I. Conventional and modified effective range parameters for the two-Yukawa model. 
The pion Yukawa was cut off at p = 500 MeV to calculate the modified ER function. In parentheses 
the ranges are expressed in terms of corresponding powers of the expected underlying scale. 



as (fm) 

ro (fm) 

rijh? (fm^) 

r 2 /A^ (fm®) 

ER 

-25.1 

1.63 (1.2/m^) 

-3.4 {-l.2/ml) 

16. (2.9/m®) 

Mod. ER 

-1.72 

0.667 (2.6/mp) 

-0.024 {-lA/mf) 

0.0028 (2.6/m®) 


The first few ER parameters from p cot 5 in Eq. (Q) are given in Table |I[ As anticipated, 
the ranges are all on the order of 1/m.,^. The modihed ER expansion leads to values more 
naturally expressed in powers of l/mp. Once again we hnd that the modified ER parame¬ 
ters truly have the long-distance physics removed, allowing for a clean £t to our EFT and 
resulting in an improved radius of convergence. Introducing some error into the value of Ott, 
as studied in the last subsection, leads not only to an inability to £t three constants but 
also to unnatural modihed ER parameters. This effect may provide a way to work towards 
a more accurate value for 

Using the conventional ER expansion to £t one or two constants in the EFT, though, 
gives results just as good as from using the modihed ER expansion. It only fails at three 
constants. Since is weak, the importance of the long-distance potential, as 

well as the need for the modihed ER expansion, are reduced. The same result occurs in the 
model of the exponential well if the strength of the potential is decreased. This behavior will 
be explained in the next section. Regardless, the robustness of the ht to the ER parameters 
for the two approaches is not qualitatively the same. The success of the conventional ER 
ht depends strongly on the initial guess for the EFT constants and the weighting of the 
ht in Eq. (pT]) by the theoretical error. These hindrances, which do not play a major role 
in the modihed ER procedure, could make it difficult to ht to real data with experimental 
uncertainties. 

By increasing we should see the radius of convergence from using the conventional 
ER ht decrease while the modihed ER result is unchanged. This is illustrated by the dashed 
lines in the left plot of Fig. which show results for the conventional ER ht after enhancing 
aTT by a factor of three. The radius of convergence is only about 300 MeV for the optimal 
Ac = 400 MeV. The results from using the modihed ER procedure are insensitive to the 
strength of the long-range coupling constant. 

In summary, for a basic two-scale problem, the modihed ER procedure provides a system¬ 
atic removal of the long-distance physics, which ensures a proper, robust ht to the low-energy 
constants, even with a strong long-range potential. This often results in a better radius of 
convergence for our EFT, although the improvement is not signihcant for a weak long-range 
potential. In addition, the modihed ER parameters rU from Eq. ([Tl|) are all on the order of 
the range of the short-distance potential, l/rup. 

III. ARE NONPERTURBATIVE PIONS NEEDED? 

We now turn to the PDS regularization scheme introduced in Ref. 0, which used the 
language of Feynman diagrams. The nonrelativistic Lagrangian of Eq. (^) gives the pion- 
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nucleon vertex 


q, a 

I 

A 

I 

I 


U 

2 /, 


-T (j-q, 


(15) 


and so, after projection onto the ^S'o partial wave, one-pion exchange contributes (g = p — p') 


9a 


^ ml , / 

~1 + Inl 1 H-— 


(16) 


to the full amplitude 


An - 1 _ 47 r 1 

M 2ip M p cot 6 — ip 


(17) 


Notice that the amplitude develops a branch point in momentum at p = —imjr/2 when the 
pion Eq. (PIB] ) is included, as expected from the effective range expansion. 

The loop integrals are evaluated in dimensional regularization with the PDS scheme, 
which subtracts all poles coming from four or fewer dimensions [^. For example, the one- 
loop bubble diagram of two nucleons with external momentum p = yME becomes 


./p Y fd q i i _ M . 

V 2 / J {2n)^ E/2 + qo — /2M + ie E/2 — qo — /2M + ie An ^ ^ 

An arbitrary mass scale p is introduced to ensure that coupling constants multiplying the 
integral Eq. (p! 8 [) retain their dimension when D is taken different from four. Physical 
quantities are independent of /i, which has two consequences: (i) the coupling constants C 2 n 
in Eq. ([I3| ) must depend on p to cancel the p-dependence of the loop integrals; and (ii) p 
must be of the same order as p or smaller to ensure that the integral Eq. ([T8|) is properly 
counted in a momentum expansion. 


A. PDS Renormalization Group 


The PDS power counting requires the constants in the effective Lagrangian scale like 
1/p and 6*2 ~ 1/p^ when considering momenta p~/iS>l/as 0. This is particularly 
relevant for S'-wave NN scattering, where the scattering length is large. Furthermore, one- 
pion exchange scales like 0{pP) as in Eq. ([T 6 |) and nucleon loops scale like 0{p) as in Eq. (|T^. 
Therefore, every n-loop diagram with insertions of the contact term Cq scales like 0{p~^). 
To have a consistent power counting, they must be summed 0 , producing the leading-order 
amplitude depicted in Fig. ^ 


^-1 = 


1 + Cq^{p + ip) 


(19) 


Matching A-i to the leading-order effective range expansion then determines the p depen¬ 
dence of Co, 
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FIG. 4. Leading order diagrams in the PDS power counting with Cq at each vertex. 


FIG. 5. All |U-dependent diagrams that contribute to the /3-function for Cq to order 


C'o(/i) 


dvr 

M 


+ l/a^ 


+ 0 (/i° 


( 20 ) 


confirming the scaling Cq l//i. 

At subleading order O(p^), perturbative corrections from one-pion exchange and C 2 
insertions begin to contribute. For non-zero pion mass, a term of the form 1)2 also 
enters. Since to this order only the combination Cq + rn^D 2 appears, we absorb D 2 into the 
dehnition of Co for convenience. Matching the amplitude to this order in the effective range 
expansion 


p cot 6 = ip + 


dvr 


dTrMo 

MAlt ^ 


( 21 ) 


determines the /i-dependence of the low-energy constants. Using the explicit expression for 
the subleading amplitude Aq in Ref. we obtain 


Coifi) 

C2{p) 


dvr /1 -b 2{p - 
M V 
MCI 
Stt 


NN 




-Pn + 1/d 

glM ' 
24vr/2 


+ 


3-8i^ + 6^ 

mi 


Ibvr 


Mm^rCo V m 


3^ 


2 + 


967r^ 


{Mm^CaY 




( 22 ) 


(23) 


where 


= A + (24) 

^TT Ann 

and Ann = ~ 300 MeV. At /i = m^, these expressions reproduce the matching 

condition quoted in Ref. 0]. The expression for Cq also coincides with the leading order 
result Eq. (p0|) for this choice of p. 

Solving the renormalization group equation for Cq is another way to generate its p depen¬ 
dence. A calculation order-by-order in the number of loops is not equivalent to solving it in 
powers of momentum due to the nonperturbative nature of the power counting. Expanding 
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FIG. 6. The left plot shows the exact phase shift in degrees for the two-Yukawa model and 
the results from one (dashed) and two (dot-dashed) constants in PDS with jjL = m-,^. The right side 
shows the error in pcot 6 for these cases (solid) as well as for the cutoff regularization (dashed) of 
Fig. g 


the amplitude, which is an observable, in powers of momentum requires that the sum of all 
graphs to a given order is /i-independent. At sub leading order, the /r-dependent graphs are 
shown in Fig. g, and applying to their explicit expressions from Ref. |g] results in the 
renormalization group equation 


p- 


dCo Mfi 


dp 


Att 


C^n + 


+ -^Cn 






1 -k 2(/i - m^)/ANN 


+ ) 


(25) 


whose solution is the expression Eq. (|^) , already obtained from matching. The solution is 
consistent with the power counting for all /r, behaving like Go ~ l//i. This is in contrast to 
a renormalization group analysis based on a loop expansion, whose solution is inconsistent 
with the power counting for p > Ann [0]- Although the above discussion was just for the ^S'o 
channel, the results also hold for the ^Si channel since the pion contribution is identical in 
the two cases to this order. 


B. PDS and the Modified Effective Range Expansion 

We can now apply the relations Eq. (^)-(|^) to generate the leading and subleading 
behavior of PDS. The error plots for the two-Yukawa model from Sec. II.C using PDS are 
compared in Fig. g with the cutoff results. Matching the constants to data ensures p- 
independent observables only to the order they were fixed; the remaining error still depends 
on the choice of p. This is different from the case with a short-range potential only , since 
pions contain /^-dependent terms to all orders in the momentum expansion. However, the 
results are relatively stable for p < 450 MeV. The radius of convergence (from extrapolating 
the slopes at small p) is about p ~ m^, which is much smaller than from cutoff regulariza- 
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FIG. 7. A Feynman-diagram identity for the Lippmann-Schwinger equation in dimensional 
regularization. 

tion.|^ In fact, using a purely short-range effective potential with cutoff regularization gives 
almost the same results as PDS with pions. 

This implies the pion is not being properly taken into account, as can be seen in the 
power counting of C 2 . Without pions, the scaling is C 2 ~ but Equation shows 

the pion contributions begin to dominate for fi ~ m^r, leading to 

C2{p > m^) -^ . (26) 

The power counting for C 2 P^ therefore breaks down for momenta above mir. A nonper- 
turbative treatment of pions might remedy this situation. It is instructive, therefore, to 
investigate how the PDS power counting differs from the modihed ER expansion, which was 
so successful in the previous section. 

First, we show that, for dimensional regularization, solving the Schrodinger equation with 
the potential Eq. ( 0 ) is the same as using the modihed ER expansion. This can be shown 
diagrammatically for the truncated potential Vsip, q) = Cq + C 2 {p^ + 5 '^)/ 2 . The amplitude 
A is found by solving the Lippmann-Schwinger equation, which is equivalent to solving 
the Schrodinger equation. In dimensional regularization, factoring the potential out of the 
loop integrals produces an additional contribution to the potential SVs{p,p) = —a-„mT^MC 2 , 
leading to Rs = V 5 -f 5Vs- Denoting the long-range pion Yukawa^ {Vl) with a double dashed 
line and the short-range delta-function terms (V 5 ) by a dot, the amplitude A can then be 
shown to satisfy the Feynman-diagram identity in Fig. ^ 

The hrst group of terms on the right-hand side of Fig. is just the one-pion exchange 
amplitude Al, which in turn is related to the long-range phase shift 5l by an expression 


^ There is a narrow window around p = 400 MeV in which the radius of convergence increases 
to about 300 MeV. Fine-tuning p gives results qualitatively similar to fine-tuning D 2 as a third 
parameter as done in Ref. |^. In that fit, however, the accuracy in the determination of the ER 
parameters and rg was sacrificed in order to minimize global errors for p < 200 MeV. In contrast, 
Eqs. (|2^ -(23) can be used to fine-tune p while still observing the matching requirements, leading 
to proper error plots and systematic results. We will not exploit this effect but instead take p = 
as in Ref. [iTtl. 


^Note from Eq. (14) that part of the pion contribution is proportional to a delta-function and is 
therefore absorbed in Vs- 
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analogous to Eq. ([T7|) . The coefficient of Vs in the denominator of the last term is the 
Green’s function (see Appendix B) 


GEir = 0, r' = 0) 


( 0 | 


E - He+ ie 


| 0 ) 


M fL{p,r) 

47rr El{p) 


(27) 


and the coefficient of Vs in the numerator is given by the square of [§] 

(28) 

Therefore the diagrammatic identity can be rewritten as 

47r |^l(p)P ^ ^Vs{p,p) ^29) 

M pcot{6 - 6 l) - ip 1 -Vs{p,p)Ge{0,0)' 

Calculating G£;(0, 0) in the PDS scheme with nonperturbative pions (see Appendix B) leads 
to the modified ER function Eq. ([TT|), 

=-—777 - - - p-a^M\n—- (30) 

MVs{p,p) ^TT 

1 1~ 2 

= + • • • , 


where and therefore are both regularization scheme dependent. Matching the hrst two 
terms in the expansion gives 


Goip) 

C2{p) 

V 


dvr 

2nro 


27rro 

CX’j^Tri'jY 


/i + In — - ( -Cjr 

T/lyj- \ Qjg 


(31) 

(32) 


where Go includes the short-distance pion contribution [see Eq. (0)]. These expressions are 
very different from the perturbative pion results Eq. (^^^) and lead to 


G2{p) ~ 


27rro 1 


(33) 


which does not exhibit a breakdown in the scaling. 

Removing the effect of the pions nonperturbatively, as in the above modihed ER proce¬ 
dure, leaves only a residual short-range potential. PDS reproduces the modihed ER expan¬ 
sion by construction in this case, with //-independent results just as for a purely short-range 
potential [^. The error plots for the two-Yukawa model, shown in Fig. |^, illustrate the 
improvement of the PDS radius of convergence from about 150 MeV with perturbative pions 
to about 800 MeV with nonperturbative pions. 

Taking pions into account perturbatively, as in the original PDS power counting scheme, 
leads to a breakdown in the scaling of the EFT constants. This is rehected by the radius of 
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FIG. 8. The left plot shows the exact phase shift in degrees for the two-Yukawa model and 
the results from one (dashed) and two (dot-dashed) constants in PDS with nonperturbative pions. 
The right side shows the error in p cot 5 for these cases (solid) as well as for the perturbative pion 
case (dashed) of Fig. 

convergence still set by as in Fig. Implementing nonpertnrbative pions, either throngh 
using the modihed ER expansion or solving the Schrodinger equation, rectihes the situation, 
leading to an improved radius of convergence set by the next scale in the underlying theory. 
It is not clear if some modihcation of the perturbative pion counting scheme can be devised 
to transcend this problem. For example, a global htting procedure that does not enforce 
an exact match of the low energy data could lead to a better overall error in Fig. |^. 
However, the ability to determine a radius of convergence graphically is lost in this case, 
requiring further analysis to determine the breakdown scale. 




C. Cutoff Regularization and the Modified Effective Range Expansion 

We return to cutoff regularization within the context of the above arguments. In this 
case, the potential Vs cannot be factored out of the loop integrals in Fig. |^, so the Schrodinger 
equation and modihed ER expansion are not as simply related. In particular, the hnite cutoff 
generates additional terms on the right-hand-side of Fig. ^ that contribute to the contact 
interactions. For example, 

) 0 (- O - + • • • . (34) 


which can be interpreted as a shift of the short-range potential Vs in the modihed ER 
expansion. The combination in parentheses originates from Eq. ( 0 ). It was shown in 
Ref. that an unnatural cutoh leads to an unnaturally large Cq, reducing the radius of 
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TABLE II. Conventional and modified effective range parameters for S-wave np scattering. 
The pion Yukawa was cut off at p = mp/2. All numbers are in appropriate powers of fm. In 
parentheses the ranges are expressed in terms of possible underlying scales, with mu = 500 MeV. 




as (fm) 

ro (fm) 

ri/h? (fm^) 

r 2 /A^ (fm^) 

% 

ER 

Mod. ER 

-23.4 

-1.70 

2.63 (1.9/m^) 
3.10 (7.9/m^) 

-0.25 (-0.091/m3) 
1.5 (24./m3) 

1.7 (0.31/m5) 

-0.53 {-bh./ml) 


ER 

Mod. ER 

5.39 

-7.46 

1.75 (1.2/m^) 
2.12 (5.4/m^) 

0.21 (0.074/m3) 

0.32 (5.2/m3) 

0.16 (0.030/m5) 

0.28 (29./m^) 


convergence of the EFT. Thus an analogous effect should follow if is large. This is why 
a strong long-range potential such as in the exponential and delta-shell model of Sec. II.A 
has better results in cutoff regularization from using the modihed ER expansion as opposed 
to the Schrodinger equation. However, a weak coupling, as occurs in NN scattering, leads 
to almost no difference in the radius of convergence. 


IV. APPLICATION TO DATA 

Now that we have a better understanding of what sets the breakdown scale for the EFT, 
we can apply the techniques of the previous sections to the problem of NN scattering. We 
take our data from the Nijmegen partial wave analysis |^, since it gives the most stable 
results over a large momentum range. As we saw in the two-Yukawa model of Sec. II.C, the 
ER and modihed ER parameters rehect the corresponding scales of the underlying physics 
quite nicely. Table || makes the same comparison using the np scattering S'-waves.Q Note 
that the mass scale needed to make the modihed ranges of order unity is actually smaller 
than the ranges r„, in stark contrast to the clean systematics of the two-Yukawa model. 

This difference can also be seen in the error plots for the ^S'o partial wave in Fig. p. As 
anticipated by the two-Yukawa model, the modihed ER and conventional ER results for 
cutoff regularization give identical results, due to the weak coupling of the pion. However, 
the optimal cutoff is Ac = 250 MeV with a resulting radius of convergence of about 300 MeV, 
consistent with the analysis of Ref. . This is much lower than the 500-600 MeV radius of 
convergence of the two-Yukawa model. The PDS results for p = (left plot of Fig. ^ break 
down earlier than cutoff regularization with a worse absolute error, although p can again 
be tuned to improve the PDS error plots (as in Sec. H.C). Using PDS with nonperturbative 
pions improves the radius of convergence (right plot of Fig. leading to results comparable 
with the cutoff, as also seen in the last section. 

The fact that the modified effective ranges and the radius of convergence reffect a smaller 
scale than 500 MeV must be due to a feature of the data that our two-Yukawa model did 
not simulate. Retardation and relativistic effects are not included in the model, but at these 
energies they are small and should be negligible. In the np system, electromagnetic effects 


^The precise values in Table || are somewhat sensitive to our choice of data and fitting procedure. 
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FIG. 9. The ^S'o np error in pcoi5 for cutoff regularization with Ac = 250 MeV (solid) and 
for PDS (dashed), both shown for one and two constants. The PDS results are for p = m-,^ and 
perturbative pions in the left plot and for nonperturbative pions in the right plot. 


are very small in phenomenological models that parameterize the data Details of 

the hadronic structure are only expected to appear for momenta above 500 MeV [§. Errors 
in the value of or in the data have some effect on the values of the higher-order range 
parameters, but do not appreciably change the ranges quoted in Table |1| or the error plot 
of Fig. |. 

There are two features which are not as easy to rule out: the pion-production inelastic 
threshold, at about p ~ 360 MeV, and correlated two-pion exchange, which gains strength 
above p ~ 400 MeV, both of which are close to the breakdown scale in Fig. Conventionally, 
the latter effect has been modeled using a a meson with a mass of about 500 MeV |^. Naive 
dimensional analysis shows that two-pion effects are the same order as the C 2 contact term, 
and could prove to be important @]. This would imply a breakdown at about m„l2 ~ 
250 MeV and modihed effective ranges with behavior r„ ~ l/mo-. The error plots of the 
data (Fig. P) are consistent with this effect, but the ranges, as shown in Table |1|, are not. 

We can test the dependence of our analysis on two-pion contributions by replacing the 
single short-distance Yukawa in the two-Yukawa model by two short-distance Yukawas rep¬ 
resenting an attractive a meson of mass 500 MeV and a repulsive vector meson (normally 
associated with both the uj and p mesons, but we will generically label it with a p)\ 


p-mpT p-rriar p-m-^r 

V(r)^^V2^Q,--_ (35) 

^ rp p ^ ^ 

Our two-Yukawa model Eq. (|T^ averaged the two short-distance effects into a single (at¬ 
tractive) short-ranged Yukawa for simplicity. Choosing an arbitrary value for a^, we tune ap 
to obtain the large scattering length of NN scattering S'-waves. We consider two representa¬ 
tive solutions: one with weak couplings ap) = (1.5,1.5) and one with strong couplings 
(ocr, ttp) = (7,14.65). The strong couplings are similar to couplings used in the Bonn poten- 


18 















tial to describe NN scattering |T^ and are the size expected from naive dimensional analysis 

M- 

When the weak conplings are used, the modihed ER approach provides a clean removal 
of the short-distance physics in both the error plots and the r„’s, similar to what was 
seen for the two-Yukawa model in Sec. II.C. Applying the same procedure to the strong 
coupling solution, however, gives remarkably similar results to those found with the data. 
In particular, the effective ranges behave similarly to those in Table y, with the underlying 
scale for the modihed effective ranges seemingly worse than the conventional effective ranges, 
and the radius of convergence drops from 600 MeV to 300 MeV in the error plots. In our 
model, we can explicitly prove that the breakdown is associated with (and not 

for example) by increasing the scalar mass. As m^, increases to rup, the results smoothly 
shift to those of the two-Yukawa model. 

Therefore, the pion is successfully removed by the modihed ER expansion, leaving the 
next scale mo-/2 to set the radius of convergence. Just as a weak pion did not require the 
modihed ER expansion for an improved radius of convergence, a weak scalar coupling a„ 
does not ahect the breakdown. This explains why only the strong coupling solution fails 
at ma/2. A strong coupling also means the range of the interaction is not 0{1) x l/rria 
I^ . The dimensionless strength of the potential is rj„ = ~ 13. Investigation of 

the three-Yukawa model Eq. (|35|) indicates that the modihed ehective ranges should scale 
as Tn ~ The results in Table || are indeed consistent with this scaling. 

Thus the large modihed ehective ranges in NN scattering imply that two-pion contri¬ 
butions are strong. If this is so, these two-pion ehects should be included in the modihed ER 
procedure as long-distance physics in order to go beyond p ~ 300 MeV. It is straightforward 
to test this conclusion in the three-Yukawa model by including the a in the long-distance 
potential. As expected, this leads to an improved radius of convergence. The next step 
should be to use the explicit form for the two-pion physics in an EFT analysis of data with 
the modihed ER procedure. 


V. SUMMARY 

An ehective held theory (EFT) approach to low-energy QCD systematically treats the 
physics of diherent length scales. For applications to nonrelativistic NN scattering and nu¬ 
clear structure, distance scales corresponding to momenta up to several times the pion mass 
may be relevant. In this regime, the long-distance physics is dictated by the pion through 
chiral symmetry. The short-distance physics is replaced by a generic expansion in powers 
of the momentum over a scale determined by the remaining underlying interactions. This 
conceptually simple procedure promises a systematic description of two-nucleon processes. 

The nonperturbative nature of nuclear forces complicates this EFT approach, leading 
to open questions about including pions within various regularization and power counting 
schemes. In this paper, we have examined cutoff regularization and dimensional regulariza¬ 
tion with the PDS scheme. Toy models were used to analyze both approaches in a controlled 
fashion before confronting NN scattering data. We hnd that the implementation of these 
schemes can be subtle. 

Within cutoff regularization, simply adding a one-pion-exchange potential to a short- 
range effective potential does not necessarily improve the radius of convergence of the EFT. 
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This is because the cutoff mixes long-distance contributions with the short-distance con¬ 
stants, which can lead to unnaturally large corrections. Only if the long-distance coupling is 
weak, as occurs for one-pion exchange in NN scattering, are these extra contributions negli¬ 
gible. By using the modified effective range (ER) procedure, however, a consistent removal 
of the long-distance physics can be implemented for all strengths of the interaction. 

In the PDS scheme, the long-distance pion physics is accounted for perturbatively with 
a consistent power counting prescription. A systematic error analysis of a two-scale model 
showed, however, that an EFT with the PDS power counting breaks down earlier than cutoff 
regularization. This breakdown is manifested by the renormalization-group scaling of the 
short-distance constants. The modihed ER procedure was used to include one-pion exchange 
nonperturbatively in dimensional regularization with PDS. This removes the long-distance 
contamination from the short-distance physics and produces error plots comparable to cutoff 
regularization. It is not clear whether a consistent power counting with perturbative pions 
can be developed to also include the benefits of this nonperturbative approach. 

With the understanding of what causes the EFT breakdown in model problems, we ap¬ 
plied the modihed ER techniques to the analysis of NN scattering data. Despite the clean 
removal of one-pion exchange effects, we found a radius of convergence of only p ~ 300 MeV 
and unexpectedly large modihed ehective ranges. We were only able to reproduce this char¬ 
acteristic behavior using a more realistic model of nuclear forces from one-boson-exchange 
phenomenology. This model includes strong attractive and repulsive Yukawa potentials re¬ 
sulting from the exchange of scalar and vector mesons. The mass of the scalar meson, which 
is associated with correlated two-pion exchange, sets the breakdown scale of the EFT, and 
the strong coupling is the source of the large modihed ehective ranges. This implies that 
to extend the radius of convergence beyond 300 MeV, we must include correlated two-pion 
ehects explicitly in the EFT. Since these ehects are strong, a modihed ER analysis could be 
essential. 
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APPENDIX A: NUMERICAL TECHNIQUES 


The two quantities needed for numerical computation of the modihed ER are the dost 
function J^{p) and the derivative with respect to r of the dost solution at the origin f'{p,r). 
Both of these can be found by solving for the logarithmic derivative of f{p,r), 


^ 1 df{p,r) 

f{p,r) dr 


(Al) 


By diherentiating Lp{r) once and using the fact that f{p,r) is a solution to the Schrodinger 
equation, we arrive at the diherential equation 
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dLp{r) 

dr 


Lp{oo) = ip, 


(A2) 


= 2mred^(r) - / 


Ll{r), 


which can be solved as a coupled equation for both the real and imaginary parts of Lp{r) 
and for Lp{0) = Lp in particular. The real part of Lp is needed in the modified ER function, 
Eq. O, and the imaginary part can be shown to be 

= (A3) 

which can be solved for the magnitude of the dost function. 

The phase of the dost function is the negative of the phase shift for the potential V{r). 
This can be found conveniently and accurately using the variable phase method 

dYr) = — sin^ [pr + d(r)] , d(0) = 0 , (A4) 

p 

with the full phase shift given by d = d(cxo). When numerically solving for a potential that 
includes a delta-shell, 

Vs{r) = -^5{r-R) , (A5) 

2??7.red 

this method must be modified to take into account the discontinuity of the wavefunction 
derivative at the radius R, 


4>'{R+) - 4>'{R-) = -A0(R_) . (A6) 

This can be accomplished by integrating Eq. ( |A4| ) out from the origin to r = R_ and then 
using the new value 




—pR + arctan 


t{R^) 

1 - 


to continue the integration to inhnity. 


t{r) = tan (pr -|- d(r)) 


(A7) 


APPENDIX B: GREEN’S FUNCTION IN THE PDS SCHEME 


We follow Ref. in solving the S-wave radial Schrodinger equation. Since the short- 
range potential Eq. (|T^ is proportional to a delta-function, we can separate it out to give 

{Hl - E)i;{r) = -Vs{r)i;{r) = -Vs{r)i;{0 ). (Bl) 


Away from the 
solutions of the 


origin, the full wavefunction -0 only depends on the regular and irregular 
purely long-range potential 


'ipiir) 


M-P) 

2ipr 


fL{p,r) + c.c. 


1 

2JEi(p)r 


fL{p,r) + c.c. 


B 'ipR{r) , 


(B2) 
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with B a momentum-dependent coefficient that is inconsequential to our result. These 
wavefunctions satisfy {Hi — E)'ipR = 0 and {H^ — E)'ipi = {Att/M) 5^{r). The Green’s 
function GE{r, r' = 0) = (0|(i? — + ie)~^\r) can therefore be written as a combination of 

the regular and irregular solution that cancels the incoming spherical wave, 


GE(r,0) 


M 

dvr 




ip \ 

Mp)\V 




M fLip,r) 
dvrr 


(B3) 


with the last equality following from Eq. (|B2|) . 

The divergent quantity G£;(0,0) can be separated into a momentum dependent term 
needed for the modihed ER function Eq. (|TT|) and a regularization scheme dependent con¬ 
stant 


- -gReGS'(0,0) = Re 


/l(p.O) 

Mp) 


reg 


—Reor(0,0) ■ 


(B4) 


The first term is the logarithmic derivative of the dost solution /l(p, r) at the origin and, 
up to a constant, does not depend on the regularization scheme. Using the PDS scheme, 
the second term is 

ReG™S(0,0) =p + a^Mln-^ + c^ . (B5) 

M 

Since a momentum-independent term can be shuffled between the two terms on the right- 
hand side of Eq. (p^, the constant is arbitrary. The imaginary part of Ge( 0, 0) can be 
read directly from Eq. (|B3|) , leading to the full Green’s function 

reg 

+ p + a-j^M In ——h Ctt + I ^ \ \2 ’ (^6) 

'fR-K \'>L\P}\ 

as used in the text. 


dvr 


- jgGJ“"(0.0) = Re 


/i,(P.O) 

reip) 
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